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We introduce a generalized angular spectrum representation for quantized light beams. By using our formalism, 
we are able to derive simple expressions for the electromagnetic vector potential operator in the case of: a) 
time-independent paraxial fields, b) time-dependent paraxial fields, and c) non-paraxial fields. For the first 
case, the well known paraxial results are fully recovered. © 2008 Optical Society of America 
OCIS codes: 000.1600, 270.0270. 



Propagation of nonclassical states of the electromag- 
netic field is an issue of growing interest in quantum op- 
tics for both fundamental and technological purposes. 1 
Consider, for instance, the relevance of propagation of 
entangled photons to quantum cryptographic systems. 2 

The purpose of the present Letter is to provide a 
novel perfectly general formalism for the representation 
of quantized light beams which can be used in any regime 
of propagation. Such objective is achieved by using a 
dispersion relation that some of us have recently in- 
troduced in Ref., 3 and a generalized angular spectrum 
representation for the field operators. 4 The usefulness 
of our approach 5 becomes manifest whenever one deals 
with quantum systems for which both paraxial and non- 
paraxial regimes of propagation may be relevant as, e.g., 
down-converted photon pairs. 6 

Consider the plane-wave expansion of the positive- 
frequency part of the electromagnetic vector potential 
operator A(r, t) = A' + '(r, t) + H.c. in the Coulomb 



gauge' 



A( + '(r,t) = J d 3 k 



167r 3 e c|k| 



1/2 



e (A) (k)o A (k) cxp (ik ■ r - ic[k|t). 



A=l 

(1) 

Since we want to describe fields propagating mainly 
along the z axis, we find it convenient to define k z = s£, 
where s = sign(/c z ) = ±1, and £ > 0. Then we can 
rewrite Eq. 0) as 

1/2 



A(+)(r,t) = ^|d^dfc^°°dC (j^ 



r3 e n c|k s 

J2 e( A )(k s )a A (k s ) exp (ik s • r - ic\k s \t), 



A=l 

(2) 

where we have defined k s = {k x ,k v , sQ). The field an- 
nihilation and creation operators satisfy the canonical 
commutation relations 7 

[a A (k s ),4,(k' s ,)] =S X y5 ss ,5( 2 \ q -q')S(t-C). (3) 
Let to > be an arbitrary frequency; at a later point 



in this Letter we shall identify to with the carrier fre- 
quency of a paraxial field. We perform a change of vari- 
ables {k x , k y , C} — * {q x , Qy,^>} such that 



k i 



C = /(q,<o), 



(4) 



where /(q, u>) > is an (almost) arbitrary function to 
be determined. For reasons that will be soon clear, we 
require /(q, ui) to increase monotonically for increasing 
u in the domain 



I u (/,q) = { W £l+:/(q, W )>0} 



(5) 



This condition implies that d/(q, u)/dw > for cu € 
X w (/, q). Therefore, in such domain, we can write 

5^ ( q _ q ') S (C - CO = 5 (2) (q - q')« [/(q, - /(q «/)] 



<5 (2) (q-q') 



6(v - to') 
d/(q,w)/dw' 



If we substitute Eq. |JBJ into Eq. © , we obtain 

[a A (k s ),4,(k' s ,)]-<5AA'^^ (2) (q-q') 



(0) 



d/(q,w)/dw 



(7) 

Equation (7J suggests the introduction of the "angular- 
spectrum" field operators a As (q, ui) defined as 



a Xs (q,uj) ee a x (k s ) \/d/(q, Ljj/duj. (8) 
By using Eq. JSJ and Eq. Q it is easy to see that 

[a Xs (q, u;), fit , s , (q', «/)] = 6\\/ 5 SS , 6^ (q - q')d(w - w'). 

(9) 

Equation © is the first main result of this Letter; it 
worthwhile to note that it is exact. No approximations 
were made to obtain it. 

The condition /(q, ui) > defines a volume V qjaJ (/) 
in the half-space R 2 x R + spanned by {q x , q y , ui]. 
This volume is bounded by the surface Sq tU (f) = 
dVq :LU (f) defined by the equation /(q, ui) = 0. 
If we define the 2-dimensional domain C q (/, ui) = 
{{lx,q y ) G R 2 : /(q,w) > 0} then, from Eq. © it read- 
ily follows 



d 2 q / da; = 



du / d 2 q. (10) 



1 



We use this equality to rewrite Eq. immediately in 
the new variables {q x ,q y ,uj} as 

1/2 



■e c 



A=l 



exp ( iq • x + isfz — \tc\J q 2 + f 2 



(11) 

where q 2 = q 2 x + g 2 , f = /(q,w), and x = Note 
that the square root of the Jacobian J = df/duj was 
used to pass from the original operators ctA(k s ) to the 
angular-spectrum operators a\ s (q, ui) . Since we want to 
develop a formalism suitable for both non-paraxial and 
paraxial light beams, we rewrite Eq. (jl 1|> as 

/>oo 

A( + )(r,i)= V / due-^ l - t - sz '^i> s {Y,t-oj), (12) 
s =±i Ja 

so that uj determines the plane carrier wave, and 
^ s (r,t]Lo) is the envelope field which, at this stage, is 
not required to be spatially and temporally slowly vary- 
ing: 



d 2 q 



fid//dw 



1/2 



X 
A=l 

x exp 
x exp 

Moreover, since k s = 



167T 3 £ <Vg 2 + / 2 < 

X] e i A) (q^)aA S (q,w) 



CM,") 

2 



iq • x + izs(f — u>/c)] 
-it (c^q 2 + f 2 -LU 



(13) 



k s \/g 2 + / 2 , where k s 



(gq + sfi,)/ yjq 2 + f 2 , and q = q/g; we have defined 



ei X \q,Lj) = e (A) (q, sf), where e { s z> (q, ui) = si x q, and 



(2), 



e W(q, W ) = (/q- S gz)/ V / g J T7 2 . 



(14) 



Until now, we furnished expressions for the field opera- 
tors in the angular spectrum representation, but not for 
the energy, the momentum, etc. However, closed expres- 
sions for these physical quantities can be easily found by 
noting that the product 

a\{k s )ax(k s ) dC = a\ s (c[, w)a As (q, ui) dui, (15) 

is invariant with respect to the change of variables Eq. 
Q. Then, for example, starting from the well known 
expression for the Hamiltonian operator of the electro- 
magnetic field (see, e.g., Ref. 7 ), after a straightforward 
calculation one obtain 

fl = \H T duJ I d 2 qft C \/g 2 + / 2 (q,c) 
1 s =±i Jo J cMM 

2 

x [ a L(q^)"A S (q,w) + a As (q,w)4 s (q,w) 



A=l 



(16) 



Similar calculations can be easily done for the other 
relevant quantities. Equation H16(l shows that, as ex- 
pected for an arbitrary field, the frequency to of the 
carrier plane wave is not equal to the frequency c|k| = 
c-J q 2 + / 2 (q, ui) of the plane- wave mode exp(ik-r). How- 
ever, as we shall see later, c|k| reduces to ui in the parax- 
ial limit. 

At this point the function /(q, ui) is still undeter- 
mined, therefore we can exploit this freedom by imposing 
some constraints on the envelope field ^ s (r, t; ui) which 
is, until now, perfectly general. In particular, we want 
to find an expression for the envelope field in which the 
Fresnel propagator 6 plays a role even beyond the parax- 
ial regime. To this end, we proceed as in Ref. 3 and we 
require *S? s (r,t = Q;u>) = 4? s (r;uj) to satisfy the time- 
independent paraxial equation: 

d 2 * s (r-Lo) , d 2 * s (r;ui) , ui d* s (r;ui) 

2is = 0. (17) 



dx 2 



dy 2 



dz 



In this way we obtain an expression for A<+)(r,f) 
which is an exact solution of the full d'Alembert equa- 
tion for any time t > and its corresponding enve- 
lope field * s (r, t = 0;uj) satisfies the time-independent 
paraxial wave equation at t = 0, as initial condi- 
tion. If we substitute from Eq. i|13|) the plane wave 
exp [iq ■ x + izs(f — ui/c)} into Eq. i|17|) . wc easily find 



/(q,w) 



g 2 c 2 
2ui 2 



(18) 



For q £ Cq(f, ui) this function satisfies all our require- 
ments: it is positive and d/(q, ui)/dui = (l + ?9 2 ) jc > 
where we have defined 8 d = qc/ (V2ui) . It is easy to see 
that the plane-wave frequency c|k| = ui(l + i? 4 ) 1 / 2 re- 
duces to ui in the paraxial limit d <C 1. Finally, a closed 
expression for the field operator A*- + ' (r, t) can be given: 



/>oc 

A( +) (r,i) = V / duje 

. ui JO 



-\Lj(t—SZ J c) 



±1 

d\( 



h{\ + d 2 ) 



'cm,") 

2 

A=l 



x exp 



V 167r 3 eocw-\/T 
( s x) (q,uj)axs((l,uj) exp 

[-iwi (yi 1) 



1/2 



i? 4 



q 2 c 

iq • x — is z 

2u) 



(19) 

Equation (|19|) is the second main result of this Let- 
ter. It is easy to recognize in the exponential func- 
tion in the third row, the sought Fresnel propagator 
in momentum space. The spatial behavior of the enve- 
lope field is entirely governed by this term. It worth 
to note that Eq. H19|l is exact, that is it has been 
obtained without any approximation and, therefore, it 
holds for both non-paraxial ($ < 1) and paraxial (d -C 
1) beams. In the latter case, the slowly varying term 
exp \—iu)t (yi + i? 4 — l)] shows that the envelope field 



2 



\l/ s (r, t; w) cannot be strictly monochromatic for any 
t > 0. 

In the remaining part of this Letter, we give two differ- 
ent examples of the application of our theory in order to 
illustrate its generality. As a first example, let us gener- 
alize the previous case and require ^ s (r,t;u>) to satisfy 
the time- dependent paraxial wave equation, for any t: 9 



d 2 * s 

dx 2 



d 2 * & 
dy 2 



c oz 



2i- 



dt 



= 0, 



(20) 



where i$! s = 1 5 r s (r, t\u>) for short. If we substitute from 
Eq. (|13|) the relevant term exp[iq • x + \zs(f — ui/c)] x 
it(c\J q 2 + f 2 — us)] into Eq. I|2(J|I . we obtain a new 



exp|— it(c-\J q* + f 2 — lo 
dispersion relation 



f(q,u) = - 1- 



2 2 

q c 

Auo 2 



(21) 



Once again, for q € C q (/, w) this function satisfies 
all our requirements: it is positive and d/(q, u))/duj = 
(l + rj 2 ) /c > 0, where rj = qc/{2uo). It is easy to see that 
the plane-wave frequency c|k| becomes c|k| = io(l +?/ 2 ). 
Also for this case a closed expression for the field oper- 
ator A( + ) (r, t) can be given: 



poo 

A«(r,i) = V / dco 

s=±l J 



d 2 q 



2 



-ia;(t — sz/c) 



1/2 



167r 3 £ocw 

X! e i A) fa' w )^ As fa' w ) ex p 



A=l 



iq • x 



.g 2 c, 

1 sz ■ 

4w v 



ct) 



(22) 

This expression is quite simpler than Eq. (|19|l . However, 
its exponential part (in the last row), differ by a factor of 
1/2 from the Fresnel propagator expression. Once again, 
we stress that Eq. Q22[l is exact, no approximation were 
made. 

As a last example of application of our formalism, we 
choose to determine the function /(q, to) by requiring w 
to coincide with the plane-wave frequency 4 c|k| : ui = 
c|k|. It is easy to see that in this case we have 

1/2 

/(q,w) = ~(l 



1 2 



? ? 

q c 

2^" 



(23) 



where the last approximate equality holds in the parax- 
ial limit qc/uj <C 1. For q G C q (/, ui) this function is 
positive and d/(q, lo)/cIlo = l/(c^/l — (qc/ui) 2 ) > 0, 
therefore all our requirements arc fulfilled. As expected, 
in the paraxial limit Eq. I|23[) coincides with Eq. Q18[l. 
Since by definition ( = \k z \ = (u)/c)\ cos 6\, it follows 
that | cos#| = /(q, ui)c/u>, and we can write 



A<+>(r,t)= £ / 



x 



cos ( 



\ 16it 3 Eocuj 



1/2 2 



/C q (,/» v— ■ - u — / A=1 

x exp [iq • x — iszu>(l — | cos6\)/c]. 



X} e s A) fa> w )aA S fa> w ) 



(24) 



Equation (|24|l is our last result. It gives an exact expres- 
sion for the electromagnetic potential vector operator of 
a generic, non-paraxial light beam, in the angular spec- 
trum representation. By expanding in Taylor series the 
| cos 9 1 term around 9 = 0, it is easy to see that Eq. 
<|24|) reduces to the well know classical paraxial expres- 
sion (with the quantum operators a\ s (q, lo) substituted 
by the corresponding classical amplitudes). Moreover, at 
the lowest order in 9, it coincides with Eq. (|19J) calcu- 
lated at the lowest order in 

In conclusion, in this Letter we presented a novel for- 
malism for the representation of arbitrary quantized light 
beams. First, we introduced an angular spectrum rep- 
resentation for the field annihilation and creation op- 
erators. Then, we used our formalism to derive an ex- 
act expression for the "paraxial-likc" envelope field of 
a light beam. Finally, we illustrated the generality of 
our theory, by applying it to the description of time- 
dependent, paraxial and non-paraxial, light beams. It 
worth to note that, although our formalism is fully quan- 
tum, all the previous results can be straightforwardly ex- 
tended to classical fields just by replacing the quantum 
operators a\ s (q, ui) with the corresponding classical am- 
plitudes a\ s (q,ui). 

We acknowledge support from the EU under the IST- 
ATESIT contract. This project is also supported by 
FOM. 

References 

1. M. I. Kolobov, Rev. Mod. Phys. 71, 1539 (1999). 

2. N. Gisin, G. Ribody, W. Tittel, and H.Zbinden, Rev. 
Mod. Phys. 74, 145 (2002). 

3. A. Aiello, and J. P. Woerdman, Phys. Rev. A 72, 
060101 (R) (2005). 

4. J. Visser, Operator description of the dynamics of optical 
modes, Ph.D thesis, Leiden University (2005). 

5. Very recently, the same dispersion relation we first intro- 
duced in Ref., 3 has been used by: G. F. Calvo, A. Picon, 
and E. Bagan, Phys. Rev. A 73, 013805 (2006), to il- 
lustrate some interesting properties of photon angular 
momentum. 

6. L. Mandel and E. Wolf, Optical Coherence and Quan- 
tum Optics, first edition, (Cambridge University Press, 
1995). 

7. R. Loudon, The Quantum Theory of Light, 3rd ed., (Ox- 
ford University Press, Oxford, 2000). 

8. See Ref. 3 for a full discussion about the physical meaning 
of the parameter 

9. I. H. Deutsch, and J. C. Garrison, Phys. Rev. A 43, 2498 
(1991). 



3 



